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Anomalous NMR Magnetic Shifts in CeCoIn 5 
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We report 115 In and 59 Co Nuclear Magnetic Resonance (NMR) measurements in the heavy fermion 
superconductor CeCoIns above and below T c . The hyperfine couplings of the 115 In and 59 Co 
are anisotropic and exhibit dramatic changes below 50K due to changes in the crystal field level 
populations of the Ce ions. Below T c the spin susceptibility is suppressed, indicating singlet pairing. 
PACS Numbers: 74.70.Tx, 76.60. Cq 



In heavy fermion systems the interplay of magnetism 
and superconductivity gives rise to a diverse range of 
ground states including an unconventional form of su- 
perconductivity. The recently discovered family of heavy 
fermion compounds CcMIxi5, where M = Co, Rh or Ir 
exemplifies these effects. Whereas the Rh compound un- 
dergoes a transition from antifcrromagnetic to supercon- 
ducting under pressure B, the Ir ^ and Co || com- 
pounds superconduct at ambient pressure, with the Co 
system exhibiting the highest known transition temper- 
ature (2.3K) for any heavy fermion system. Evidence 
from heat capacity, thermal transport and /iSR indicate 
that the pairing symmetry in the superconducting state 
is unconventional and that there are line nodes in the 
superconducting gap. 

The bulk magnetic susceptibility, x, of tetragonal 
CeMIns displays systematic trends consistent with the 
diversity of observed ground states. In all three cases x 
is anisotropic, and is largest for field applied along the 
c direction. In the ab plane, Xab is essentially the same 
for all three materials. However, Xc exhibits a maximum 
at - 10 K for CeRhIn 5 (T N = 3.8 K), whereas for the 
superconductors Celrlns and CeCoIns Xc diverges at low 
temperatures until T c is reached. For both of these ma- 
terials Xc also exhibits a plateau-like feature around 50 
K, which is less pronounced for the Ir system. The origin 
of this feature and the relationship between Xc and T c 
have been sources of debate, however both the plateau 
and the divergence are intrinsic and independent of field. 

I 

Here we report a detailed study of site-specific mag- 
netic shifts in CeCoIns using nuclear magnetic resonance 
(NMR). Measurements in the normal state provide a mi- 
croscopic measure of the local susceptibility and we find 
anomalous temperature dependencies. This behavior is 
likely due to the thermal depopulation of a crystal field 
(CEF) excitation of the Ce ions. We find remarkably 
strong departures from the expected proportionality be- 
tween bulk susceptibility and the NMR Knight shift. We 
will argue that this effect is indicative of a high degree 
Ce moment localization, a feature that may play a role in 
the mechanism for superconductivity in this material. In 
the superconducting state the temperature dependencies 
of the shifts reveal a suppression of the spin susceptibility 
consistent with spin-singlet pairing. 



Crystals of CeCoIns were grown from an In flux as 
described in S. The tetragonal crystal structure of 
CeCoIns consists of alternating layers of Celn3 and CoIn2 
and so has two incquivalcnt In sites per unit cell. The 
In(l) site has axial symmetry and is analogous to the sin- 
gle In site in cubic Celn3. There are four low symmetry 
In(2) sites per unit cell, two on each of the lateral faces of 
the unit cell, located a distance 0.306c above and below 
the Ce-In layer. § The zero field 115 In (I = 9/2) Nu- 
clear Quadrupolar Resonance (NQR) spectrum reveals 
an axially symmetric site with 115 i/q(1)= 8.173 ± 0.005 
MHz, and 77(1) = 0.0 at 4 K, whereas the electric field 
gradient (EFG) at the In(2) site is characterized by 
115 i/ Q (2)= 15.489±0.001 MHz, and nC2) = 0.386±0.001, 
where vq and r\ are defined as in |8]||. The NMR spec- 
trum of the 59 Co (/ = 7/2) indicates a site with axial 
symmetry and 59 vq = 234±1 kHz at 4 K. Both the In 
and the Co EFG's are essentially temperature indepen- 
dent, varying less than 0.5% between 4 K and 100 K, 
indicating that significant structural changes are absent 
in this temperature range. 

The magnetic shift measurements were made on a 
large single crystal of CeCoIns, which was mounted with 
the c axis either parallel or perpendicular to the ex- 
ternal field, for fields between 3 and 5 T. Field-swept 
spectra were obtained by measuring the spin echo in- 
tensity as a function of applied field at fixed frequency. 
The shifts were determined by measuring several of the 
115 In transition fields -ff cxp for each site at several dif- 
ferent fixed frequencies. The nuclear spin Hamiltonian 
H = (hv Q /6)[3I 2 - I 2 + V (I 2 - I 2 )} + 7 RI • (1 + K) ■ Ho, 
where K = (K a , Kb, K c ) is the magnetic shift tensor, was 
diagonalized and the resonance fields H res for each tran- 
sition and each In site were then calculated. The spectra 
were then fit by minimizing x 2 = ^2i(H re s — H oxp ) 2 as a 
function of (6, <j>, K a ,Kb, K c ), where 9 and 4> are the po- 
lar angles relating Ho to the crystal axes (a, b, c). Note 
that such a procedure is necessary because the strong 
quadrupolar interaction gives rise to a significant angular 
dependence of H res so that even a misalignment of 1-2° 
can cause a significant error (~ 30%) in K. The Co shift 
and EFG were determined by measuring the positions of 
the central and satellite transitions at fixed field. 

Given three nuclei and two possible field orientations 
for each there are seven distinct magnetic shifts. Note 



that for In(l) and Co the magnetic shift is isotropic in 
the ab plane, whereas for In(2) the shift differs for Ho 
parallel or perpendicular to the unit cell face. The tem- 
perature dependencies of K for both In sites as well as 
the Co are shown in Fig. (1), together with x f° r both di- 
rections. K is a measure of the local electronic spin den- 
sity at the nuclear site. In general, the shift is given by: 
K(T) = K + J2i AXi{T), where K is an orbital shift, 
independent of the local spin density at the nuclear site 
and the temperature, and is the hyperfine coupling to 
Xi, the i th component of the susceptibility x — YliXi- 
Both Kq and Ai can be anisotropic. All of the magnetic 
shifts except 115 K (2) j_ are proportional to x f° r T > 40K 
for Ho||a6 and T > 60K for Ho||c. Below these tempera- 
tures n5 K(2)n, 59 K ab , 59 K C and n5 K(l) c show dramatic 
departures from x- Furthermore, K(2)x is not pro- 
portional to x m an y temperature regime, and exhibits 
a dramatic downturn below 40 K [note that the axis for 
115 K(2)j_ is reversed in Fig. (1)]. Figure (2) shows K 
versus x f° r both field directions. 
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FIG. 1. The magnetic shift versus temperature. The solid 
lines show Xc and Xab- The dashed line is a fit to \ as de- 
scribed in the text; the dotted lines, which have been offset 
vertically to coincide with the shift scales, are calculations as 
described in the text. [Note the reversed axis for 115 K (2) ±.] 

Note that K oc x f° r high temperatures (T > 40 K), and 
the intercept and slope give Kq and A, whose values are 
listed in Table (1), where Aht is determined for high 
temperatures, and Alt for low temperatures. The Co 
shifts track those of the In(2) for both directions, where 
A c (Co)/A c (In(2)) = 0.26 and A ab (Co)/A ± (ln(2)) = 
0.33. Therefore, it seems likely that the Co is not di- 
rectly coupled to the Ce, but couples to the Ce only via 



the In(2). 

Anomalous departures from K oc x have been known 
to exist in Ce compounds for several years, although the 
reason for the departure is still under debate. fl0|-|l2|| It is 
generally considered that the Ce 4f electron does not have 
a significant direct overlap with the orbitals of neighbor- 
ing nuclei. Rather, it is the 6s and 5d orbitals of the Ce 
that are hybridized, and the Ce 4f moment can create a 
hyperfine field at a neighboring atom by polarizing the 
conduction electrons at the Ce site, which is then trans- 
ferred to the neighbor via an RKKY interaction. The 
conduction electrons at the neighbor then create a hyper- 
fine field at the nucleus via a contact interaction. Two 
different mechanisms have been proposed to explain the 
anomalous shift behavior in other heavy fermion systems. 
In CeSn3 K (Sn) and x differ below ~ 150 K, and this ef- 
fect has been ascribed to modifications of the effective 
hyperfine coupling at the Sn (via the RKKY interaction) 
by the onset of Kondo compensation below a tempera- 
ture Tk- ]i"l| , |i3| in CeCu2Si2 Ohama and coworkers ob- 
served that the Cu and Si magnetic shifts also exhibit 
departures from K oc x below ~ 100 K, and they at- 
tribute this behavior to the depopulation of an excited 
CEF level of the Ce ions (J = 5/2) and not Kondo co- 
herence. jl2l In this case, the overlap between the Ce 4f 
orbitals and the conduction electrons differs depending 
on the CEF level populations, resulting in temperature 
dependent hyperfine couplings to the Cu and Si. 
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FIG. 2. K versus \ in both directions. The dotted lines 
are linear fits to the high temperature data, and the symbols 
are the same as in Fig. (1). 

In fact, the measured shifts in CeCoIns show behaviors 
similar to those observed in CeCu2Si2. Namely, the K 
versus x plots exhibit positive slope at high tempera- 
tures; however, at low temperatures K oc x IS recov- 
ered, but with a negative slope. Ohama et al. attribute 
the negative hyperfine coupling to an orbital overlap be- 
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tween the ligand s orbital and the Ce 4f orbital. They 
distinguish this direct transferred hyperfine mechanism 
from that in which the 4f moment polarizes the con- 
duction band at the Ce site. According to Ohama et 
al, the direct contribution can become negative when 
only the lowest CEF doublet is occupied. Heat capac- 
ity data in CeCoIns suggest the presence of an excited 
CEF doublet at ~ 50K above the ground state doublet, 
H so it would be reasonable to ascribe the anomalous 
shift behavior in CeCoIns to depopulation of an excited 
CEF doublet. The strong site dependence of the shift 
anisotropy in CeCoIn 5 also suggests a direct coupling be- 
tween the In or Co and an anisotropic, localized Ce 4f 
orbital, as in CeCu2Si2- Recent deHaas van Alphen and 
photoemission data in CeCoIns are consistent with LDA 
calculations, which assume the Ce 4f electron is itinerant. 
flit In contrast, the strong directionality of the hyperfine 
couplings in these materials implies an electronic struc- 
ture that is more tight binding rather than free electron 
like. In both systems the temperature scale for Kondo 
compensation is much lower than the CEF splitting, sug- 
gesting that more localized 4f character, which could be 
a source for spin fluctuations, is essential for the devel- 
opment of heavy fermion superconductivity. 

In order to investigate the possible role of CEF ef- 
fects we have fit x to extract the CEF parameters. The 
dashed lines in Fig. (1) show a fit to the expression 
X 1 = Xcef + where xcef is the CEF suscepti- 
bility for the Ce ion, and A is a molecular field term. 
The Ce ion in CeCoIns experiences a crystal field with 
tetragonal symmetry, so Hcef = b^O® + + b\0\, 
where the O™ are the Steven's operators. J15J In this 
field the J = 5/2 manifold is split into three dou- 

.(1) T^2) 



Tf), where the wavefunctions are given 



biets (r 6 ,r^ 

by: ±|)> Tsina|±|) =F cosa|=F|), =Fcosa|±|) ± 
sina|=F§), and xcef = (d 2 (logZ)/dH 2 ) H=0 . Here Z 
is the partition function for the Hamiltonian Ticc = 
Wcef + S./MsH • J, where gj — 6/7, is the Bohr 
magneton, and J is the spin operator for J = 5/2. We 
find the best fit for the T 6 ground state (J z = ±^), 
with excited states at 34K and 102K above the ground 
state, a = 1.47, and an anisotropic molecular field: 
A c = 18.8 mol/emu and X a b = —113.2 mol/emu. 0] 
The anisotropy of A reflects Ce-Ce couplings which differ 
for neighbors in and out of the plane. The fit repro- 
duces the plateau feature, and suggests that the anoma- 
lous behavior of the magnetic shifts below 50K may also 
be explained by changes in the hyperfine couplings as 
the excited CEF states are depopulated. Note, for exam- 
ple, that in Fig. (2) 115 K(1) C appears to be independent 
of Xc at low temperatures. This behavior suggests that 
115 K(1) C couples only to the excited CEF states. If we 
decompose xcef = Xcef + Xcef m *° contributions from 
the ground state doublet and from the excited doublets, 
one might then expect K = Kq + A gs x es + A cx x c * , where 
(xT 1 = (Xcef)" 1 + A. We determine Xcef (Xcef) b y 
suppressing the field dependence of the excited (ground) 



state energy levels in the expression for Z. By adjusting 
K , A gs and A ox appropriately, we can qualitatively ex- 
plain the temperature dependence of all the shifts in Fig. 
(1) (dotted lines). Note that Xcef < in the ab plane, 
so for this component the absolute value of X a b was used. 

The anomalous behavior of the shifts might also be 
explained by two components of x with a different ori- 
gin than crystal field states. However, there is only one 
Ce site in the unit cell, and susceptibility and heat ca- 
pacity data indicate that the observed properties can be 
entirely attributed to the Ce (i.e., Co is nonmagnetic in 
CeCoIns). Therefore it seems likely that the two com- 
ponents can only be attributed to different CEF states 
on the Ce ions. It is interesting to note that measure- 
ments of the In(l) shift in the isostructural compound 
CeRhlns reveal a positive hyperfine coupling for 4 K 
< T < 50 K, with no signs of the dramatic departure 
from K oc x seen in CeCoIns. Jl?]] Clearly, if the hyper- 
fine anomaly is the only mechanism at work in CeCoIn 5 
then the CEF parameters in CeCoIns must differ signif- 
icantly from those in CeRhIn 5 . In fact, recent work by 
Takeuchi and coworkers suggests that the ground state 
CEF level in CeRhIn 5 is T 7 rather than r 6 . ITJ] 
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FIG. 3. The magnetic shift for Ho||ab in the superconduct- 
ing state. The solid line is the bulk susceptibility, which be- 
comes fully diamagnetic below 2.3 K. 

Below T c x is dominated by the diamagnetic response, 
masking the intrinsic behavior of the spin susceptibility. 
K, however, couples only to the spin susceptibility and 
provides a direct measure of ^ spm in the superconduct- 
ing state. The temperature dependencies of the shifts 
for both In sites as well as the Co in CeCoIns are shown 
in Fig. (3) for Ho[|a6 down to 1.4K . Because of the 
thin platelet morphology of CeCoIns, demagnetization 
fields in the superconducting state can be significant for 
H ||c, precluding an accurate determination of the mag- 
netic shift since the local field at the nucleus is poorly 
determined. We estimated that for our sample, the de- 
magnetization factor for H||c is N c /Att « 0.79. Thcrc- 
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fore, although we observe a decrease in the resonance 
frequencies for this direction one cannot resolve whether 
the decrease is due to a change in K or to a change in 
Hq internally. However, for Ho||a6 the demagnetization 
factor is much smaller, so the internal field below T c is 
known to a greater degree of accuracy. Therefore, we 
only present data on the shifts for the field in the plane. 



In(1) 




H ±c 

H = 33.3 kG 



Co 




H ±c 

H = 39.0 kG 
J i i l_ 



ln(2) 




H ±c 

H = 31.5kG 





40.0 40.1 39.72 39.80 39.6 39.8 
Freq (MHz) 

FIG. 4. Spectra of the In(l), Co and In(2) at various tem- 
peratures through T c . The In(2) spectrum is for Ho normal 
to the unit cell face. The two peak structure of the In(l) spec- 
trum indicates the presence of two slightly differently oriented 
crystals in the sample. 

The decrease in 115 K (l) a & seen in Fig. (3) implies a de- 
crease in x spin - However, 59 K ab , 115 K(2) ± and 115 K{2) 1{ 
increase below T c [note the reversed axis for 115 K(2)± in 
Fig. (3)]. Spectra of the In(l), In(2) and Co at different 
temperatures are shown Fig. (4), clearly exhibiting the 
behaviors seen in Fig. (3). An increase of the absolute 
value of 115 K (2) ±, 115 K (2) n , and 59 K ab below T c can be 
understood by recognizing that the hyperfine coupling is 
negative below 50 K [see Fig. (3)], so an increase in K 
implies a decrease in x spin . Thus, all of the shifts for 
H 1 c are consistent with a decrease in x splI \ implying 
spin-singlet pairing of the Cooper pairs in the supercon- 
ducting state. Given the recent heat capacity and ther- 
mal conductivity measurements revealing higher orbital 
symmetry, [^],[l9| we can conclude that the order parame- 
ter in CeCoIns has d-wave symmetry. During the course 
of this work, we became aware of similar work by the 
group of Kohara |2(J who report magnetic shift results 
below T c . Although our conclusions about singlet pairing 
are the same, the temperature dependencies of the shifts 
differ. 

We thank S. Dunsiger for assistance with the measure- 
ments, as well as D. MacLaughlin, R. Heffher, and J. 
Lawrence for valuable discussions. This work was per- 
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TABLE I. The hyperfine couplings and orbital shifts of the 



In(l), In(2) and Co. 


shift a 


K {%) 


Art (kOe/p B ) 


Alt (kOe//i B ) 


In(l)c 


0.79(5) 


8.94(34) 


-0.4(1) 


In(l) oi 


0.13(4) 


12.08(40) 


12.08(40) 


In(2) c 


-2.10(3) 


32.4(3) 


22.8(3) 


In(2) N 


0.76(2) 


10.26(17) 


-12(1) 


In(2) ± 


1.10(1) 





-34.7(9) 


Co c 


1.00(1) 


8.4(5) 


6.20(5) 


COai, 


0.68(1) 


3.30(9) 


-4.20(19) 
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